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LOCAL GROUP DWARF ELLIPTICAL GALAXIES: II. STELLAR KINEMATICS TO LARGE RADII IN 

NGC 147 AND NGC 185 
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ABSTRACT 

We present kinematic and metallicity profiles for the M31 dwarf elliptical (dE) satellite galaxies 
NGC 147 and NGC 185. The profiles represent the most extensive spectroscopic radial coverage for 
any dE galaxy, extending to a projected distance of eight half-light radii (8r c g ~ 14'). We achieve 
this coverage via Keck/DEIMOS multislit spectroscopic observations of 520 and 442 member red giant 
branch stars in NGC 147 and NGC 185, respectively. In contrast to previous studies, we find that both 
dEs have significant internal rotation. We measure a maximum rotational velocity of 17 ± 2 km s _1 
for NGC 147 and 15 ±5 km s _1 for NGC 185. While both rotation profiles suggest a flattening in the 
outer regions, there is no indication that we have reached the radius of maximum rotation velocity. 
The velocity dispersions decrease gently with radius with an average dispersion of 16 ± 1km s~ x for 
NGC 147 and 24± 1 km s" 1 for NGC 185. The average metallicity for NGC 147 is [Pe/H] = -1.1 ±0.1 
and for NGC 185 is [Fe/H] = —1.3 ±0.1; both dEs have internal metallicity dispersions of 0.5 dex, but 
show no evidence for a radial metallicity gradient. We construct two-integral axisymmetric dynamical 
models and find that the observed kinematical profiles cannot be explained without modest amounts 
of non-baryonic dark matter. We measure central mass-to-light ratios of M/Ly = 4.2 ± 0.6 and 
M/Ly = 4.6 ± 0.6 for NGC 147 and NGC 185, respectively. Both dE galaxies are consistent with 
being primarily flattened by their rotational motions, although some anisotropic velocity dispersion 
is needed to fully explain their observed shapes. The velocity profiles of all three Local Group dEs 
(NGC 147, NGC 185 and NGC 205) suggest that rotation is more prevalent in the dE galaxy class 
than previously assumed, but is often manifest only at several times the effective radius. Since all 
dEs outside the Local Group have been probed to only inside the effective radius, this opens the 
door for formation mechanisms in which dEs are transformed or stripped versions of gas-rich rotating 
progenitor galaxies. 

Subject headings: galaxies: dwarf — galaxies: kinematics and dynamics — galaxies: individual 



(NGC 147, NGC 185) 



1. INTRODUCTION 

Dwarf elliptical (dE) galaxies are characterized by low 
surface brightness (/U e ff,v > 22 mag arcsec -2 ), little to 
no gas, and old to intermediate age stellar populations. 
These galaxies are spatially clustered and account for 
more than 75% of objects brighter than My < —14 
in nearby galaxy clusters ([Ferguson fc Sandagel HMD : 
iFer guson fc Binggelil[T99l . Since there are few, if any, 
isolated dE galaxies, most proposed dE galaxy forma- 
tion scenarios are environmentally-driven, such as via 
gas stripping or gravitational harass ment of gas-rich spi- 
ral or dwarf irregu l ar galaxies (e.g. [Dekel fc Silk) 119861 : 
iMoore et all 119981 : lAguerri fc Gonzalez-Garcial (2009( 1. 
Recent evidence for embedded stella r disks and bars in 
some dEs suppo r t these scenarios dJerien et all |2000|; 
Geha et al.l 120051: iLisker et all 120061: iChilingarian et all 



2007f ). Since gas-rich systems have significant rotation 



velocities, a key test of this hypothesis class is compar- 
ing the kinematics of dE and gas-rich galaxies. 
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Internal kinematic measurements for dEs in the For- 
nax and Virgo Cluster have revealed an unexpected di- 
chotomy: roughly half of the observed dEs have sig- 
nificant rotation about the major-axis, while the re- 
main ing galaxies show no detectable major-axis rota- 
tion (iPedraz et al.ll2q02t iGeha et al.l[200l Ivan Zee et "all 
120041 : Ide Riicke et al.ll2005( l. dE galaxies with no rotation 
present a challenge to any scenario i n which dE progen - 
itors are gas-rich and thus rotating (Mo ore et al.l ll998). 
Ram pressure stripping alone does not affect the kine- 
matics of the collisionless stellar component. Gravita- 
tional processes such as galaxy harassment can reduce 
the amount of rotational angular momentum, converting 
this into random velocity dispersion, however, simula- 
tions are unable to completely remove rotatio nal support 
(|Maver et al.ll2001t iMastropietro et al.ll2005D . 

Because dEs have low surface brightnesses, integrated 
light (long-slit) kinematic observations of dEs in both 
galaxy clusters and the Local Group are currently lim- 
ited to within the half-light effective radius (r c g). Thus, 
it is possible that significant rotational angular momen- 
tum lies at larger distances than has so far been explored. 
Addressing this question requires kinematics tracers at 
larger radius. While the dynamics of dE globular clus- 
ter systems suggest that there may be signifi cant rota- 
tion at large radius (jBeaslev et al . 2006, 2009) , the small 
number of clusters limits the usefulness of this tracer. 
The dEs in the Local Group share the same general 
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properties as the more numerous dEs in galaxy clusters, 
howe ver, their proximit y allows us to resolve individual 
stars ([Geha et al J 12006). and thus trace kinematics pro- 
files down to arbitrarily low surface brightness and large 
radius. 

The Local Group contains three dE galaxies, NGC 205, 
NGC 147 and NGC 185, which in addition to the com- 
pact elliptical M32, are the brightest satellites around 
the spiral galax y M31. We fo l low th e terminology con- 
ventions of, e.g. iBender et al.l (JT992), although alterna- 
tives are in use. While the differing terminologies can 
be confusing, the key point is that NGC 205, NGC 147, 
and NGC 185 on the one hand, and M32 on the other 
hand, should not be grouped in a single class. The former 
three galaxies have fundamental plane properties that in- 
dicate similarity to galaxies in the dwarf spheroidal class 
(dSph). By contrast, galaxies such as M32 have proper- 
ties that indicate similarity to (giant) elliptical galaxies. 
This likely indicates a differe nce in formation history , as 
discussed more fully in, e.g.. iKormendv et al.l (|2009l) . 

The prototype of the dE galaxy class is NGC 205 
which is in very close projection to M31 (40' = 8kpc). 
Both photometric and kinematic evidence su ggests that 
NGC 205 is tidal l y int eracting with M31 (jChoi et al.l 
120021: IGeha et al.l 12006ft and may in fact be on its 
first orbital approach (jHowlev et al.l I2008D . In con- 
trast, NGC 147 and NGC 185 lie at a projected dis- 
tance of 7° (~ 150 kpc) from M31. While the metal- 
poor stellar halo of M31 extends to this distance 
(jGuhathakurta et alJl2005t iKalirai et al.H2006f >. its gravi- 
tational influence of M 31 is not expected to influence the 
kinematics of these dE satellites (see ? I4.5[) . NGC 147 and 
NGC 185 are separated from each other by merely 58', 
leadi ng to speculation that they may be a bound galaxy 
pair (jBaadd ri944t Ivan" den Bcrgh 1998). An accurate 
distance estimate suggests that the se two galaxies are 
physi cally separated by over 60 kpc ()McConnachie et all 
I2005D and may not be bound (however, see § 14. 5p . The 
large observed distances between these two dEs and their 
parent M31, combined with their relative proximity to 
the Sun, make NGC 147 and NGC 185 the best available 
targets for studying the properties of dE galaxies to large 
radius. 

NGC 147 and NGC 185 share several similar funda- 
mental properties, such as absolute luminosity (My ~ 
— 15.5), half-light radius (r e g ~ 2' = 0.3kpc), and central 
velocity dispersion (for exact values see Table 1). How- 
ever, they differ markedly in many aspects. NGC 185 
contains some gas, dust and evi dence for recent star for- 
matio n confined to its center (Martmcz-Dclgad o et al.l 
[l999h . while NGC 147 is devoid of gas or dust and shows 
no si gn of recen t star formation activity (jYoung fc Lol 
H997t lSage et alJll998f >. An aver age metallicity of [Fe/Bfl 
= —1.4 is inferred for NGC 185 (jMartmez-Delgado et all 
H99l . while a more metal ric h population of [ Fe/H] = 
-1.0 is inferred for NGC 14 7 (IHan et al.l 119971: IDavidgd 
12001 IGoncalves et~alll2007h . HST/WFPC2 imaging of 
both galaxies also implies the presence of intermediate- 
age stars, with NGC 147 ha ving a more significant con- 
tribu tion than NGC 185 (Bu tler fc Martmez-Delgadol 
2005). Thus, while both dEs are dominated by old to 
intermediate-age stars, the metallicity and age mixture 
of these components are different in each galaxy. 

The first spectroscopic long-slit observations of Local 



Group dEs by IBender et al.l (|1991l) suggested that while 
NGC 185 appeared supported entirely by random mo- 
tions, NGC 147 had a significant rotational component 
of 6.5 ± 1.1km s - 1 . Som ewhat deeper observations by 
iSimien fe Prugniel ()2002h extended the kinematic pro- 
files out ~ 2' (lr c ff), but did not measure a significant 
rotation for either galaxy. In contrast, both studies mea- 
sured a rotational component in NGC 205 (within the 
tidal radius) , later confirmed by IGeha et ail (|2006l ) to 
be 11 ± 5 km s" 1 . U sing the Simien & Prugniel data, 
Ide Riicke et al.l (|2006l ) constructed dynamical models for 
all three Local Group dEs, concluding that, within the 
half-light radii, all three galaxies have mass-to-light ra- 
tios above that expected from the stellar populations 
alone. 

In this paper, we present accurate radial velocities 
out to 8r e g in the Local Group dE galaxies NGC 147 
and NGC 185, based on Keck/DEIMOS multi-slit spec- 
troscopy for 520 and 442 red giant branch (RGB) stars, 
respectively. The paper is organized as follows: in §[2]we 
discuss target selection for our DEIMOS slitmasks, the 
observing procedure and data reduction. In §[31 we dis- 
cuss the velocity, velocity dispersion and metallicity pro- 
files. In §|4l we construct two-integral dynamical models 
for these galaxies and discuss the best fit mass distribu- 
tion and dark matter content. Finally, in §[5j we discuss 
the implications of these results for the formation of dE 
galaxies. 

Througho ut this paper we ad o pt di stance moduli de- 
termined bv lMcConnachie et al.l ([2005D via the tip of the 
RGB method. The distance modulus for NGC 147 is 
(m - M) = 24.43 ± 0.04 (675 ± 27 kpc) and NGC 185 
is (m - M) = 24.23 ± 0.03 (616 ± 26 kpc). This places 
NGC 147 and NGC 185 at projected distances of 140 and 
185 kpc from their parent galaxy M31, respectively. 

2. DATA 
2.1. Target Selection 

Individual stars in NGC 147 and NGC 185 were tar- 
geted for spectroscopy according to their probability of 
being a RGB star based on Canada-France-Hawaii Tele- 
scope CFH12K mosaic ima ging in the R- and /-bands 
(jBattinelli fe Demersll2lolallH ). The CCD mosaic covers 
a 42' x 28' region centered on each galaxy with 0.206" 
pixels. Stellar photometry from these images extends 
2 mag below the tip of the RGB. The spatial distribution 
of stars is shown in Figure [T] Color- magnitude diagrams 
are shown in Figured] 

The absolute magnitude of the RGB for a metal-poor 
stellar population is Mi ~ —4. At the assumed distance 
and reddening of NGC 147 and NGC 185, the apparent 
magnitude of the RGB population is I a ~ 20.3 (see Fig- 
ure [2|) . Since internal metallicity variations will cause 
a spread in the colors of RGB stars, we select spectro- 
scopic targets based primarily on apparent magnitude 
rather than color. We note there is a 0.1-0.2 mag pho- 
tometric zero-point error between individual CCD chips 
in the CFHT photometry which we discovered after our 
spectroscopic observations were completed. While this 
led us to use an alternative dataset to determine the sur- 
face brightness profiles in § 14.11 the error had little to no 
affect on our targeting efficiencies. 

For both dEs, highest priority in the spectroscopic tar- 
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Fig. 1. — Palomar Sky Survey images of the two Local Group dE galaxies, NGC 147 (left) and NGC 185 (right). Images are 30' X 30'; 
North is up, East is to the left. The major-axis of each galaxy is indicated by the solid line. The spatial position of Keck/DEIMOS 
spectroscopically confirmed member stars are shown as blue circles, non-members are plotted as open squares. 



get list was assigned to stars between 20.5 < I < 21.0. 
Second priority was given to stars between 20.0 < I < 
20.5 and 21.0 < I < 21.5; lowest priority was assigned to 
objects I > 21.5. To minimize Galactic foreground con- 
tamination, targets were required to have {R — I) Q > 0.2. 
Stars with photometric errors larger than 07 > 0.1 
were rejected, as were stars having neighbors of equal 
or greater brightness within a radius r < 4". The photo- 
metric criteria used for spectroscopic target selection are 
shown in Figure [2] 

Slitmasks were created using the DEIMOS dsimulator 
slitmask design software. For each slitmask the software 
fills in the mask area to the extent possible with the high- 
est priority input targets. It then fills in the remaining 
space on the slitmask with lower priority targets. The 
specifics of each mask are listed in Table 2. 

2.2. Observations and Data Reduction 

The data were taken with th e Keck II 10-m tel escope 
and the DEIMOS spectrograph (jFaber et al.ll2003f >. Four 
multislit masks were observed in NGC 147 on the nights 
of August 28 and September 5-7, 2005. Five masks were 
observed in NGC 185 on September 16, 2006. Mask po- 
sitions, exposure times and other observing details are 
given in Table 2. The masks were observed with the 
1200 line mm -1 grating covering a wavelength region 
6400 — 9 100 A. The spectral dispersion of this setup is 
0.33A, and the resulting spectral resolution, taking into 
account the anamorphic distortion, is 1.37A (FWHM, 
equivalent to 47km s" 1 at the Ca II triplet). The spatial 
scale is 0.12" per pixel and slitlets were 0.7" wide. The 
minimum slit length was 4" which allows adequate sky 
subtraction; the minimum spatial separation between slit 
ends was 0.4" (three pixels). 

Spectra were reduced using a modified version of the 
spec2d software pipeline (version 1.1.4) developed by the 
DEEP2 team at the University of California-Berkeley 
for that survey. A detailed descripti on of the two- 
dimen sional reductions can be found in Simon fc Gehal 
(20071). The final one-dimensional spectra are rebinned 
into logarithmic wavelength bins with 15 km s _1 per 



pixel. An example DEIMOS spectr um taken with th e 
same settings is shown in Figure 4 of Ge ha et ail (|2006f) . 

2.3. Radial Velocities and Error Estimates 

Radial velocities were measured by cross-correlating 
the observed science spectra with a series of high signal- 
to-noise stellar templates. Stellar templates were ob- 
served with Keck/DEIMOS using the same setup as de- 
scribed in § !2.2l and cover a wide range of stellar types (F8 
to M8 giants, subgiants and dwarf stars) and metallici- 
ties ([Fe/H] = -2.12 to +0.11). We calculate and apply 
a telluric correction to each science spectrum by cross 
correlating a hot stellar template wit h the night sky ab - 
sorption lines following the method in S ohn et al.l {2007). 
The telluric correction accounts for the velocity error due 
to mis-centering the star within the 0.7" slit caused by 
small mask rotations or astrometric errors. We apply 
both a telluric and heliocentric correction to all veloci- 
ties presented in this paper. 

We determine the random component of our velocity 
errors using a Monte Carlo bootstrap method. Noise is 
added to each pixel in the one-dimensional science spec- 
trum, we then recalculate the velocity and telluric cor- 
rection for 1000 noise realizations. Error bars are defined 
as the square root of the variance in the recovered mean 
velocity in the Monte Carlo simulations. The system- 
atic contribution to th e velocity error was determined by 
iSimon fc~G cha (2007j) to be 2.2 km s _1 based on repeated 
independent measurements of individual stars. The sys- 
tematic error contribution is expected to be constant 
as the spectrograph setup and velocity cross-correlation 
routines are identical. We add the random and system- 
atic errors in quadrature to arrive at the final velocity 
error for each science measurement. Radial velocities 
were successfully measured for 1227 of the 1483 extracted 
spectra across the nine observed DEIMOS masks. The 
majority of spectra for which we could not measure a 
redshift did not have sufficient signal-to-noise. The fit- 
ted velocities were visually inspected to ensure reliability. 

2.4. Membership Criteria 
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Fig. 2.— Color-magnitude diagrams for NGC 147 (left) and NGC 185 (right) based on IBattinelli fc Demeri H2004al fbT) R- and /-band 
CHFT12K photometry, corrected for extinction, in a 42' X 28' region centered on each galaxy. Blue circles are Keck/DEIMOS spectroscop- 
ically confirmed members of each galaxy, open squares are non-members. The numbered regions indicate the location of our primary (1), 
secondary (2), and tertiary (3) spectroscopic priorities. 



The expected velocity distributions of NGC 147 and 
NGC 185 overlap with both the stellar velocity distribu- 
tion of the M31 halo and the foreground Milky Way disk 
and halo. Thus, a simple velocity cut is not sufficient to 
define membership in these dE galaxies. We determine 
membership in NGC 147 and NGC 185 using a modi- 
fied v ersion of the probability method of iGilbert et al.l 
(2006). This method computes a membership probabil- 
ity for each star based on the assumption that dE stars 
have different physical properties than the contaminating 
foregr ound and background populations. 

The IGilbert et al.l (|2006l ) method computes the prob- 
ability that a given star is either a RGB star at the 
distance of NGC 147 or NGC 185, or a foreground 
Milky Way dwarf star based on the star's location in 
four spectroscopic and photometric diagnostics: line of 
sight velocity, strength of the surface-gravity sensitive 
Na II absorption line at ~ 8190A, position in a color- 
magnitude diagram, and a comparison of photometric 
(CMD-based) and spectroscopic (based on the Ca II 
triplet at ~ 8500A) metallicity estimates. A star's po- 
sition in a diagnostic is compared to probability distri- 
bution functions based on appropriate training sets to 
determine the likelihood the star is a RGB or foreground 
dwarf star. The likelihoods from the individual diagnos- 
tics are then combined to produce an overall likelihood. 

For the present analysis, sever al modifications wer e 
made to the method published in IGilbert et al.l ([2006), 
which was optimized to identify red giant branch stars 
in the halo of M31. Several of the diagnostics utilize 
Johnson- Cousins V and / magnitudes, while the dE pho- 
tometry was taken in the R and / bands. The CMD 
diagnostic and spectroscopic versus photometric metal- 
licity diagnostic were re-derived using IVandenBerg et al.l 
(2006) isochrones in (/, R — I) shifted to the distance 
of each dE. To utilize the Na II diagnostic, Vanden- 
Berg isochrones were used to define a relation between 
R — I and V — I, which was then used to derive a.V — I 
color for each star. The velocity distribution of each dE 
is modeled as a Gaussian, with a constant velocity dis- 



persion and a central velocity that varies with distance 
along the semi-major axis. We expect little contamina- 
tion f rom M31 halo stars at th ese radii (< 1 star per 
mask: iGuhathakurta et al1l2005D . 

Stars that are more than three times more likely to be 
a red giant branch star at the distance of the dE than a 
foreground dwarf star, and which have a spatial position 
and line-of-sight velocity consistent with that of the dE, 
are classified as members. Given these membership crite- 
ria, we identify 520 and 442 RGB star as members of the 
dE galaxies NGC 147 and NGC 185, respectively. Ve- 
locities and velocities errors for individual member stars 
are listed in Tables 6 and 7 in the electronic version of 
the paper, or on request from the authors. 

3. RESULTS 

The measured velocities of individual stars allow us 
to probe the dynamics of NGC 147 and NGC 185 to 
much larger radius than possible via integrated-light 
spectroscopy. In order to compare to previous results, 
and to construct dynamical models, we bin the individ- 
ual velocity measurements based on their major-axis dis- 
tance x, using bin widths Ax = 2.05' for NGC 147 and 
Ax = 1.87' for NGC 185. The bin widths were chosen 
to roughly match the effective radius of each dE. The 
bin with integer index number i spans values of x from 
(i — j)Ax to (i + \)Ax. The value of x for each star was 
calculated using spherical trigonometry and a zenithal 
projection, using the galaxy center positions and major 
axis position angles (Table 1 and 3); positive radii lie on 
the North-East side of each galaxy. 

We analyze our kinematic data following the procedure 
descri bed in Appendix A of Ivan der Marel fc Anderson! 
( 2009) . For each bin number i we calculated the weighted 
mean velocity Vi of the individual stellar velocities Vj, 
and its formal error AVi. Individual weights were set 
equal to the quadrature sum of the observational veloc- 
ity errors Avj and the velocity dispersion Oi for the given 
bin. The latter was determined by finding the value that 
maximizes the likelihood of the set of velocity residuals 
Vj — Vi. The formal error A<Ji was calculated using a 
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Fig. 3. — (Top) Individual velocities as a function of semi-major axis distance for dE member stars (blue circles) and non-members (open 
squares) of NGC 147 (left) and NGC 185 (right). Positive radii correspond to the North-East side of each dE. (Top middle) The binned 
major-axis velocity profile showing significant rotation velocities in each dE galaxy. (Bottom middle) The velocity dispersion profiles are 
nearly flat to the last measured radial point. (Bottom) The [Fe/H] spectroscopic metallicity as a function of semi-major axis distance. 
Barred errors indicate the error on the mean metallicity, non-barred errors indicate the metallicity dispersion in each bin. 



Monte-Carlo approach, taking into account the observa- 
tional velocity errors Avj . The systemic radial velocity, 
V sys , for each galaxy was determined by finding the value 
that provides the lowest x 2 when the dynamical models 
described in §2] were fit to the data and is subtracted 
from each of the binned major-axis velocities Vi . The sys- 
temic velocity for NGC 147 is V sys = -193.1 ±0. 8km s" 1 
and V" sys = -203.8±l.lkm s" 1 for NGC 185. This is con - 
sistent with literature values from [Bender et al.l ()1991[ ) 
and the velocity compilation from the NASA Extragalac- 
tic Database. 

To further improve the signal-to-noise ratio we folded 
the profiles above by taking the weighted average of bins 
on opposite sides of the galaxy (i.e., bins i and — i), using 
the expected anti-symmetry and symmetry of the rota- 
tion and dispersion profiles, respectively. Values of the 
observed V ms and their formal errors were obtained by 
adding the results for V and a in quadrature. The full 
kinematic profiles are shown in Figure |3l the folded pro- 
files are shown in Figure El 

3.1. The Velocity Profiles and Rotation in dE Galaxies 

We present the major-axis velocity profiles for our two 
dEs in Figure [3] The profiles extend to slightly more 
than 14' in each galaxy, equivalent to 7 r c g and 9 r e s in 



NGC 147 and NGC 185, respectively. Previous kinematic 
obse rvations for each dE have been limited to within one 
r er r (|Simien fc PrugnieU l2002| ) . In Figure El a clear ro- 
tation signal is seen even in the individual unbinned ve- 
locity distribution of stars. The binned major-axis ve- 
locity profiles in this figure show a well-defined rotation 
curve profile. The maximum observed rotational veloc- 
ity is 17 ± 2 km s" 1 and 15 ± 5 km s" 1 for NGC 147 and 
NGC 185, determined from the maximum data point in 
the folded rotation curves (Figure EJ ■ While our velocity 
profiles flatten at large radius, it is unlikely that, even 
at these large radial distances, we have reached the max- 
imum rotation velocity in these galaxies. The velocity 
dispersion profiles in both dEs decline gently with radius. 
The weighted average values are a = 16 ± 1 km s _1 and 
24 ± lkm s" 1 for NGC 147 and NGC 185, respectively. 

The significant rotation velocities we measure are in 
fact consistent with literature results of little to no rota- 
tion. Previous observations of NGC 147 and NGC 185 
were limited to within 2', corresponding to the in- 
ner three data points of Figure El Fitting a simple 
straight line to our inner data points, we predict a 
rotation velocity of merely 0.8 km s" 1 at this radius 
for NGC 185, we l l belo w t he observational errors o f 
both lBender et al.l (119911) and lSimien fc Prugniell (2002), 
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and fully consistent with their conclusion of no rota- 
tion. For NGC 147, we predict a rotation velocity of 
4.3 km s" 1 within 2'. This i s mar ginally inconsistent 
with both ISimien fc Prugniel (12002|) who fo und no ro- 
tation at this radius. and lBender et al.l (|1991| ) who mea- 
sured 6.5 ± 1.1 km s _1 , but is consistent with the average 
of these two observations. Our discovery of rotation at 
radii larger than r c g has significant implications for stud- 
ies of more distant dE galaxies. We will discuss this issue 
further below and in §[5j 

Early kinematic observations of the L ocal Group dEs 
(|Davies et al.l 119831: iBender et aLlll991[) established the 
paradigm that dE galaxies are supported by anisotropic 
velocity dispersions, in contrast to low-luminosity ellip- 
tical galaxies which are rotationally-supported. Our ob- 
servations show that the Local Group dEs do indeed have 
significant rotational support. To demonstrate this, we 
compare the ratio of the observed and predicted val- 
ues of the maximum rotational velocity to the average 
velocity dispersion using the quantity (V max /a)*. For 
a rotationally-supported oblate galaxy, (Vmax/c)* = 1; 
values less than unity imply anisotropic-support. Us- 
ing the measured ellipticity of each dE (Table 1), the 
measured maximum rotation velocity and the aver- 
age velocity dispersion values from above, we compute 

(fmaxAO* = °- 95 for NGC 147 and (KW<T)* = 0.91 

for NGC 185. Because our observed maximum rotational 
velocities are lower limits on the true maximum rotation, 
these (V max /<r)* values are also lower limits to the true 
values. For the third dE in the Local Group NGC 205, 
we measured (V; nax /cr)* = 0.3 in lGeha et aTl(|2006l ). The 
original value for NGC 205 may have been somewhat 
larger given that the outer parts of this galaxy have 
been tidal disrupted by interactions with M31. Thus, 
while the majority of kinematic support for the observed 
shapes of the Local Group dE galaxies comes from ro- 
tation, some anisotropic velocity dispersions are needed 
to fully explain the observed shapes. This is consistent 
with our dynamical modeling results in § 14 .31 

3.2. Metallicity 

We measure the spectroscopic metallicities ([Fe/H]) of 
individual RGB stars based on the Ca II triplet absorp- 
tion lines near A = 8500 A . This method relies on an em- 
pirical calibration to Galactic globular clust ers and is re- 
liable for metallicities above [Fe/H]> —2.5 (|Kirbv et al.l 
2008), which is the case for our member stars. We calcu- 
late the equivalent widths (EWs) of the three Ca II ab- 
sorption lines u s ing the line and continuum definitions of 
iRutledge et al.l (|1997af l and convert this to [Fe/H] using 
the em pirical calibration relationship of Rutlcdg e et al.l 
(|1997bl) . We determine the error on the [Fe/H] values 
with the Monte Carlo method described in § 12.31 and add 
in quadr ature a systemat i c unce rtainty of 0.3 A as deter- 
mined in lSimon fe Gehal (|2007|) . We restrict this analy- 
sis to stellar spectra with signal-to-noise ratio greater or 
equal to 5 per pixel. 

The binned metallicity profiles for NGC 147 and 
NGC 185 are shown in the bottom panels of Fig- 
ure [3] For NGC 147, we find an average metallic- 
ity of [Fe/H]= -1.1 ± 0.1 and for NGC 185 we find 
[Fe/H]= —1.3 ± 0.1. This is in agreement with pre- 
vious measurements from a variety of sources, includ- 



ing photometric and spectroscopic metallicity estimates 
(THan et al.lll997tlMartmez-Delgado et al.lll999l: iDavidgd 
2005 HGoncalves et al.ll2007D . While the mean metallicity 
is well determined, there is significant internal metallicity 
dispersion of 0.5 dex in each galaxy, assuming a Gaussian 
metallicity distribution function. Although the metallic- 
ity distribution function is asymmetric, given our errors 
the assumption of a Gaussian distribution is justified. In 
Figure [3l we plot both the error on the mean (barred 
error bars) in each radial bin and the internal metallic- 
ity dispersion (unbarred errors). Despite our significant 
radial coverage, we do not find evidence for radial metal- 
licity gradients. 

4. DYNAMICAL MODELING 

We have constructed a set of dynamical models to pro- 
vide an initial int erpretation of the da t a. Th e software 
was developed by Ivan der Marel et all (|1994l ) for mod- 
eling of th e Local Group galaxy M32 a n d used more 
recently in Ivan der Marel fc van Dokkuml (|2007a| ). The 
models assume a two-integral form of the distribution 
function and are used to fit the binned profiles in Fig- 
ure [5j Binning results some information l oss and meth- 
ods do exist to model unbin ned data (e.g. IStrigari et al.l 
[200l IChaname et~aT1 [200l iCoccato et al.1 12009T ). ho~ 
ever we chose to bin the kinematic data here primar- 
ily to facilitate comparisons to previous work on these 
galaxies. A full analysis of this dataset should preferably 
avoid spatial binning of the data, restrictive use of only 
the lowest-order velocity moments and allow for a more 
general form of the distribution function. We intend to 
construct such models in a fut ure paper using the d is- 
crete Schwarzschild method of IChaname et a l. (2008). 

Our present models are based on solutions of the Jeans 
equations of hydrostatic equilibrium with the following 
assumptions: (1) an oblate axisymmetric geometry with 
a constant axial ratio; (2) a constant mass-to-light ra- 
tio M/L; and (3) a two-integral distribution function 
of the form / = f(E,L z ), where E is the energy and 
L z the angular momentum around the symmetry axis. 

Two-integral models have = v%, so they can be 
thought of as the axisymmetric generalization of spher- 
ical isotropic models. A summary of the equations on 
which the models are based is provid ed in Section 2 of 
Ivan der Marel fc van Dokkuml (|2007aD . We do not re- 
peat this here, but focus instead on the specific appli- 
cation to NGC 147 and NGC 185. Quantities for these 
galaxies that are relevant to the models are listed in Ta- 
ble [3] and are discussed below. 

4.1. Brightness and Density Profiles 

The modeling starts with the observed surface pho- 
tometry of each galaxy. While the CFHT data discussed 
in § 12.11 are appropriate for measuring surface brightness 
profiles, an unrecoverable calibration error resulted in a 
significant zero-point error between CC Ds. Instead, we 
use the data presented by iKentl (|1987j ) in the Gunn r- 
band. We transformed the Kent major-axis data to the 
Johnson U-band using the (B — V) colors from the RC3 
(see Table [3]) combined with the tra nsformation equa - 
tion V - r = 0.486(B - V) - 0.273 porgensen|[l99l . 
The resulti ng surface brightne ss profiles, corrected for 
extinction (|Schlegel et al.l 119981) . are shown in Figure [4] 



NGC 147 



20 



i 22 - 



26 



I 1 1 1 I 1 1 1 I 



Local Group dE Dynamics 

NGC 185 

18 | — i — i i i i i i — r 



i ... i . i i i 



2 4 6 

x (arcmin) 




-1 -0.5 0.5 1 
log(x[arcmin]) 



4 22 - 



24 



26 



I ... I 



2 4 6 8 

x (arcmin) 




-1 -0.5 0.5 1 
log(x[arcmin]) 



Fig. 4. — Extincti on-corrected V-band surface brightness profiles in mag/arcsec 2 of NGC 147 (left) and NGC 185 (right). Black dots 
show the data from Kent (1987), transformed to the V-band. Data inside of 1.5 times the seeing FWHM were ignored. The top panels 
have the linear major-axis radius x along the abscissa, while the bottom panels have log(x) along the abscissa. Red curves in the top panels 
are the best Sersic- profile fits to the data. Red dots in the bottom panels show the Sersic extrapolation of the data to larger radii. Blue 
curves in the bottom panels are the predictions of the best-fitting luminosity density parameterization j(R, z) in equation JTJ. Parameters 
of the fits are listed in Tabic [3] 



and extend to 5.3' and 6.8' for NGC 147 and NGC 185, 
respectively. For both galaxies this is less than half the 
extent of our kinematical data. To extrapolate the pho- 
tometry to larger radii we determined the best-fitting 
Sersic profiles, which are over-plotted in the top panels of 
Figure [4] The profiles fit the data well, with RMS resid- 
uals of 0.04 and 0.08 mag for NGC 147 and NGC 185, 
respectively. 

We set the apparent axial ratio q a (q a = 1 — e) equal to 
the average of the observed axial ratios in the radial range 
1' — 5' (see Table [3]). Any systematic variations in axial 
ratio over this radial range are small in both galaxies 
(< 0.1; Kent 1987). The total apparent magnitude, Vq, 
calculated using this axial ratio and our best fitting Sersic 
model above agree with that in the RC3 catalog to within 
~ 0.1-0.2 magnitudes for each dE galaxy. This provides 
added confidence in the use of Sersic extrapolations at 
large radii. 

The inclination of each dE is not well constrained by 
the data, however, both the predicted kinematical pro- 
files and the inferred M/L are generally insensitive to the 
assum ed inclination (e.g. Ivan der Marel fe van Dokkuml 
l2007alfbT ). The distribution o f intrinsic axial ratios qt fo r 
elliptical galaxies is known (|Tremblav fc Merrittl Il995h . 



and we assume th is distribution for our dEs. We use this 
distribution as in Ivan der Marel fc van Dokkuml (2007a) 
to calculate the statistically most likely inclination. Since 
NGC 147 appears flatter on the sky than NGC 185, its 
most likely inclination is closer to edge-on (t = 78°) than 
that for NGC 185 (i = 64°). The intrinsic axial ratios 
qt corresponding to these assumed inclinations are given 
in Table [3j We discuss the dependence on inclination in 

§H31 

We parameterize the three-dimensional luminosity 
density as: 



tf + Af • (i) 



Here z is the symmetry axis of the galaxy and (R, (j>, z) 
are the usual cylindrical coordinates. We fit the pro- 
jected surface brightness (obtained by numerical line- 
of-sight integration of j) to the combination of the 
Kent (1987) data and its Sersic-profile extrapolation to 
20'. The surface brightness profile fits are shown in the 
bottom panels of Figure HI and the model parameters are 
listed in Table [3] These results are independent of the 
assumed inclination, although jo depends on inclination 
through the ratio (qt/q a ), as well as on the adopted dis- 
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tances. The surface brightness models fit the extrapo- 
lated profiles adequately, with RMS residuals out to 20' 
of 0.06 and 0.10 mag for NGC 147 and NGC 185, respec- 
tively. 

4.2. Dynamics 

Given the j(R,z) luminosity density, we assume a 
constant M/L to calculate the gravitational potential. 
We then solve the Jeans equations to determine the in- 
trinsic dynamical quantities, followed by a luminosity- 
weighted line-of-sight projection to calculate observable 
quantities. In particular, the models yield the quantity 
V Trns = [V 2 + a 2 } 1 / 2 , where V is the mean velocity and 
cr is the velocity dispersion. The quantity V lms scales 
linearly with (M/L) 1//2 , but its spatial and radial depen- 
dence are otherwise uniquely determined. The quanti- 
ties V and a can also be predicted separately, but this 
requires the introduction of a free parameter (or func- 
tion) k that specifies how the second azimuthal velocity 
moment separates into mean and random components: 

^=fcg-^] 1/2 - (2) 

For k = the model is non-rotating, whereas for \k\ = 1 
the velocity dispersion tensor is isotropic and the model 
is a so-called "oblate isotropic rotator" . The quantity k 
is therefore similar to the quantity (V/a)* discussed is 
§ 13.11 with the difference that k is defined locally while 
{V/a)* is defined in terms of globally averaged projected 
quantities. 

4.3. Data-Model Comparison 

For proper comparison to the data, the model predic- 
tions were binned on the sky in rectangles of size Aa; x Ay. 
Here y is the minor-axis direction. The quantity Ay 
was taken to be the approximate extent of the region 
around the major-axis over which data was obtained. For 
NGC 185 we also took into account that the observed 
stars are on average somewhat offset by Sy « —1.5' from 
the major-axis (Figure [T]). No luminosity weighting was 
applied in calculation of the binned model predictions, 
since the observed stars are distributed approximately 
homogeneously in each bin. 

The value of M/L was chosen to best-fit the folded 
VJ-mg profile near the galaxy center, which in practice was 
chosen to be the region inside 4'. We determine mass-to- 
light ratios of M/L v = 4.2 ± 0.6 and M/L v = 4.6 ± 0.6 
for NGC 147 and NGC 185, respectively. The value of k 
in our models was chosen to best fit the rotation curve 
at all available radii. The parameters of the best-fitting 
models are listed in Table[3J Figure[5]compares the corre- 
sponding model predictions to the observed kinematical 
profiles. We discuss the M/L values further in S I4.4I 

For NGC 147, the predicted profile provides a statisti- 
cally acceptable match to the V TmB data (x 2 = 3.8 for 8 
degrees of freedom or xl — 0.5 ). Although the M/L was 
chosen to only optimize the fit near the center, the pre- 
dicted profile fits the full observed profile out to the last 
measured data point. Hence, even data out to ~ 7.3r G g- 
provides no evidence for an increase in M/L with ra- 
dius. By contrast, for NGC 185, the predicted Kms pro- 
file does not provide a statistically acceptable match to 
the data. The fit is reasonable inside 4', since the M/L 



was chosen to optimize the fit there. However, for the 
data points outside of this radius the observed weighted 
average V rms = 23.0 ± 1.1km s _1 . The average model 
prediction at these radii is V lms = 20.6 km s -1 . There- 
fore, the observations exceed the model predictions at 
large radii with 2.2cr significance. The most natural in- 
terpretation of this result is that the M/L of NGC 185 is 
increasing with radius, although tides may also produce 
this affect (see § |4~5|) . 

The predicted rotation curves for both galaxies provide 
a statistically acceptable match to the data (x 2 = 11 -9 
and 7.6, for 8 and 7 degrees of freedom, for NGC 147 and 
NGC 185, respectively). However, inspection of the pre- 
dicted rotation curves shows that they may be slightly 
too high at small radii and slightly too low at large radii. 
This could easily be improved by choosing k to be an 
increasing function of radius. This would be no less arbi- 
trary than choosing it to have a constant value as a func- 
tion of radius, as we have done here. If the model also 
provides a statistically acceptable fit to Vr ms , as is the 
case for NGC 147, then this guarantees that the velocity 
dispersion a is fit as well. By contrast, for NGC 185 an 
increase of M/L with radius might be a more natural 
way to improve the fit to the rotation curve, given the 
V rms results. The best-fitting k values in Table[3]indicate 
that while both galaxies have significant rotation rates, 
neither galaxy rotates fast enough to account for its flat- 
tening. Hence, the velocity dispersion tensors must be 
anisotropic, in agreement with our conclusion from § 13. 11 

The fact that NGC 147 is well fit by a constant- M/L 
two-integral model does not mean that it is not embed- 
ded in a dark halo. The profile of V^ ms with radius is 
almost flat. So it is likely that models with a logarith- 
mic gravitational potential can fit the data equally well. 
Conversely, the fact that NGC 185 is not well fit by 
constant M/L two- integral model does not mean that 
it necessarily must have a dark halo. It may be possi- 
ble to fit the data with a constant- M/L model that has 
a different velocity dispersion anisotropy than assumed 
here. Alternatively, our assumption of dynamical equilib- 
rium may be incorrect for NGC 185 at large radius and 
we discuss the possibility of tidal interactions in § 14.51 
In principle it is possible to constrain the velocity dis- 
persion anisotropy independent of the assumed gravita- 
tional potential by using information on deviations of the 
line-of-sight velocity distributions f rom Gaussians (e.g., 
iCarollo et al.l Il995t iRix et aD 1 1997T ) . To further explore 
these issues, we intend to use the discrete Schwarzschild 
method of Chaname et al. (2009) to construct dynamical 
models with three-integral distribution function, includ- 
ing a dark halo (i.e., M/L varying with radius), which 
will avoid spatial binning of the data. Such models will 
not necessarily yield smaller error bars or better con- 
straints on quantities of interest (such as the mass-to- 
light ratio or presence of a dark halo), given the limited 
number of stars available to trace the kinematics. How- 
ever, they will have the benefit of yielding more rigorous 
error bars and constraints, which include, e.g., modeling 
uncertainties due to the unknown velocity anisotropy and 
radial mass-to-light ratio profiles of these systems. 

4.4. Mass-to-Light Ratios 

We compare our inferred M/L ratios with stellar 
population models to determine whether or not non- 
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Fig. 5. — Major-axis kinematics of NGC 147 (left) and NGC 185 (right). Black dots show the binned results extracted from our data as 
described in the text. Black curves show the predictions of the best-fitting f(E, L z ) models. Parameters of the fits are listed in Table l3l 
The top panels provide evidence for an M/L increase with radius in NGC 185, but not in NGC 147. 




Fig. 6. — M/ Lb versus extinction corrected color (B — V)o- Observed values for NGC 147 (red) and NGC 185 (blue) are shown as data 
points with error bars. Curves show model predictions from Chariot & Bruzual (2007) for the IMFs from Chabrier (2003; left) and Salpeter 
(1955; right). Each curve is for a fixed metallicity Z and connects the predictions of simple stellar populations of fi xed age. Solid circle s 
indicate ages of 10 Gyr for each metallicity. The predicted curves extend up to the age of the Universe (13.73 Gyr; Spcrgcl ct al. 2007). 
Black curves are, from left to right, Z = 0.0001, 0.0004, 0.004, and 0.008. The red and blue curves shows the interpolation of these curves 
to the metallicity appropriate for NGC 147 and NGC 185, respectively. We note for the preferred Chabrier IMF, the inferred M/Lg ratios 
of these dE galaxies cannot be explained via their stellar content alone. 



baryonic dark matter is needed to explain our ob- 
served kinematics. Despite our restrictive model as- 
sumptions, our resulting M/L values are likely ver y ro- 
bust. Both Ivan der Marel fc van Dokkuml (|2007bD and 
Cappellari ct al. (2006) have compared the results of sev- 
eral dynamical modeling studies, showing that different 
modeling approaches produced consistent M/L ratios to 
within a systematic accuracy of ~ 5%. This is true de- 
spite the use of varying data sets (long-slit and integral 
field data) and modeling approaches (two- and three- 



integral models, and models with and without dark ha- 
los). In the context our models, we have further explored 
the robustness of our inferred M/L values by construct- 
ing models with a range of inclinations (from edge-on to 
inclinations corresponding to models as flat as qt = 0.4). 
Consistent with previous studies we find that both the 
M/ L values and the predicted profiles of V rms with radius 
are insensitive to the assumed inclination (well within the 
error bars). We also explored models with different ways 
of extrapolating Kent's (1987) surface brightness data to 
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larger radii. Again, we found the effect on the M/L val- 
ues and the predicted profiles of Vi- ms with radius to be 
well within the error bars. 

We first com pare our M/L values to those of 
Ide Riicke etUl (|2006l) who modeled NGC 147 and 
NGC 185 using a three-integral method. For the compar- 
ison, we transform our M / Ly values to the -B-band. The 
-B-band values were obtained from the modeled U-band 
values using the galaxy (B — V)o color (Table 3) and the 
solar color (B — V)q — 0.65. The de Rijcke et al. value for 
NGC 147, M/L B = 4.0^2 4, is consistent with our result. 
However, their value for NGC 185, M/L B = 3.0±q °, is 
lower than our result at 1.5<7 significance. The underlying 
reason for this may be that NGC 185 has an increasing 
M/L with radius, as suggested by Figure [5] The best- 
fitting M/ L is then likely to depend on the radial extent 
of the data, which is larger for o ur data set than for th e 
integrated- light measurements of Ide Riicke et all (|2006[) . 

Figure [6] shows the mass-to- light ratio M/Lb for each 
dE versus the extinction corrected color (B — V)q. Wc 
compare these measurements to the stellar population 
predictions of Chariot & Bruzual ( 2007, priv. comm.) 
which is an updated version of their iBruzual fc Chariot] 
(2003) models. The [Fe/H] metallicities, determined in 
§GE3 are Z = 0.0015 ([Fe/H] = -1.1) for NGC 147 and 
Z = 0.0009 ([Fe/H] = -1.3) for NG C 185 which assumes 
a sola r composition of Z Q = 0.0189 ( [Anders fe Grevessi 
1989). The predictions at these metallicities are over- 
plotted in Figu re \§\ for b oth t he ini tial mas s func- 
tions (IMFs) of IChabrieri ((20031) and ISalpeteri ((19551) 
over the range from 0.1-100 M©. For the Chabrier 
IMF, the maximum M/Lb ratio at the mctallicity of 
NGC 147 and NGC 185 is attained f or a popula- 
tion as old as the Universe (13.73 Gyr; iSpergel et al.l 
I2007T) with values M/Ls !max = 2.5 and 2.3, respec- 
tively. This is not as high as observed, so for the 
Chabrier IMF one must invoke the presence of non- 
baryonic dark matter to fit the observed M/Lb val- 
ues. The Salpeter IMF tracks, while closer to the ob- 
servations, are only marginally consistent with obser- 
vations for very old stellar populations. Stellar popu- 
lation studies suggest signi ficant intermediate age stel- 
lar po pulations in both dEs (B utler fc Martmez-Delgadol 
2005), thus the Salpeter IMF also implies the presence 
of dark matter. The Chabrier IMF is now generally 
preferred over the Salpeter IMF, as the latter has too 
many stars at low masses to properly fit rea l populations 
(e.g.. IBell fc de Jondl2000l : IChabrierl[200l . Therefore, 
the high observed M/L values in NGC 147 and NGC 185 
must be due at least in part to the presence of dark mat- 
ter. 

Figure [7] shows Mj Lb versus the velocity dispersion 
a. For NGC 147 and NGC 185 we use along the abscissa 
the weighted average V^ ms inside 4', as listed in Table 1. 
For comparison we show in Figure [7] the sample of (gi- 
ant) elliptical galaxies compiled by van der Marel & van 
Dokkum (2007b). The latter follow a tight linear relation 
of the form M/L B oc a. The dE galaxies NGC 147 and 
NGC 185 clearly do not follow this relation. This result 
is not specific to these two dE galaxies. For comparison 
we plo t also the results for NGC 205 from iGeha et al.l 
(2006), and for six dE galaxies in Virgo from Geha et 
al. (2002). These all fall in the same region of the plot. 
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Fig. 7. — M/Lg versus velocity dispersion a. Observed values 
for NGC 147 and NGC 185 are shown in red and green respectively. 
Cyan and blue data points are other dE galaxies from de Rijcke et 
al. (2006) and Geha et al. (2002), respectively. Black data points 
are giant ellipticals from the compilation of van der Marel & van 
Dokkum (2007b). 

The same distinction between dwarf and giant ellipticals 
is also revealed by plots of the fundamental plane; see e.g. 
the plot of K3 vs. Ki in Geha et al. (2003). A traditional 
interpretation of this is that it indicates a fundamental 
distinction between dwarf and giant ellipticals. However, 
Graham & Guzman (2003) have proposed more recently 
that there is actually a continuous sequence between the 
two, with the upturn in M/L towards dE galaxies driven 
by the fact that their Sersic n parameter is lower than for 
giant ellipticals. Either way, it is clear from Figure[7]that 
dE galaxies have M/L values that are similar to those of 
intermediate- to high-luminosity giant ellipticals. The 
latter are redder and have higher metallicity (and quite 
possibly older stellar populations) than NGC 147 and 
NGC 185. A natural explanation for the observed M/L 
of NGC 147 and NGC 185 is that these galaxies have a 
significant dark matter component, consistent with our 
conclusions from the comparison to population synthesis 
models. 

4.5. A Bound Pair? 

Given our velocity and mass estimates, we revisit the 
question of whether or not NGC 147 and NGC 185 form 
a gravitationally bound pair. The angular distance be- 
tween NGC 147 and NGC 185 is 58' = lOkpc Assum- 
ing distances measured by IMcConnachie et ail ((20051) . 
NGC 147 and NGC 185 are separated by a line-of-sight 
distance of 59kpc, and thus a total physical separation 
of r = 60 kpc. For two point masses to be gravitationally 
bound, their potentia l energy must exceed the kinetic 
energy of the system (|Davis et al.l 119951 : Ivan den Bergh 
11998ft . This leads to the criterium b = (2GM sys /rAv 2 ) > 
1 for a system to be bound, where M sys is the combined 
system mass. We measure, in §[3j a radial velocity differ- 
ence between these two dEs of Av = 10.7 ± 1.3 km s _1 . 
Approximating our dE galaxies as point masses, we cal- 
culate the criteria for these two objects to be bound as 
b = 1.6 ± 1.1. This criteria suggests that, in the absence 
of transverse motion, NGC 147 and NGC 185 are a grav- 
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itation ally bound pair, in agreement with Ivan den Bergh 
41998ft . However, the large error on this quantity, in ad- 
dition to the likely non-zero transverse motions of these 
two satellites, leaves wide open the possibility that these 
are not in fact bound. 

NGC 147 and NGC 185 are satellites of M31. Our 
equilibrium models assume that the outer kinematics of 
these dEs are not affected by tidal interactions with M 31, 
or each other. This is not the case for the third dE in 
the Local Group, NGC 205, which is in the process of be- 
ing tidal stripped at radii beyond 4.5' = lkpc by M31. 
To confirm our equilibrium assumption for NGC 147 and 
NGC 185, we calculate their instantaneous tidal radius in 
the rotating frame by again approximating these galaxies 
as point sources. We as sume a total mass for M31 be- 
tween 0.8-1.6 x 10 12 M W (|Klvpin et al.ll2002HSeigar et al.1 
2008). We determine a tidal (Jacobi) radius between 10 
to 12kpc (greater than 25r e ff) in each dE. For compar- 
ison, the tidal radius of the interacting dE NGC 205 is 
0.7 kpc. The internal kinematic profiles of NGC 147 and 
NGC 185 are unlikely to be currently affected by tidal 
interactions with M31, at the distances probed by the 
data presented here. However, these calculations do not 
preclude the possibility of past tidal interactions, cither 
between these two galaxies or with M31, and it is yet 
possible that there are unbound stars in our kinematic 
sample. The regularity and bi-symmetry of our kine- 
matic data for NGC 147 suggest that any contamination 
is minimal. In NGC 185, the upturn of Vr ms profile over 
our models at large radius may either be interpreted as 
an increasing M/L ratio with radius, or possible tidal 
contamination. Deeper imaging and kinematic observa- 
tions are needed, particularly for NGC 185, to determine 
to what extent these dEs are affected by tidal interac- 
tions at radii beyond 8r c g. 

5. CONCLUSIONS 

We have presented mean velocity, velocity dispersion 
and metallicity profiles for the Local Group dE galaxies 
NGC 147 and NGC 185 based on Keck/DEIMOS spec- 
troscopic observations of 520 and 442 member RGB stars 
in each dE, respectively. The profiles represent the most 
extensive spectroscopic radial coverage for any dE galaxy, 
extending to a projected distance of eight half-light radii 
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(14' ~ 8r c ff). Contrary to previous results, we find that 
both dEs have significant rotation velocities which con- 
tributes to the observed flattened shapes of each galaxy. 
Our two-integral dynamical modeling suggests that the 
observed rotation velocities cannot fully explained the 
observed shapes and that some anisotropic velocity dis- 
persion is required. Our modeling estimates M/L ratios 
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pected from stellar populations alone. Thus, some dark 
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TABLE 1 

Local Group dEs at a Glance 



Name 


a (J2000) 
(h : m : s) 


8 (J2000) 
(»:':") 


Type 


Dist. 

(kpc) 


m v 


A v 


M vfi 


€ 


^eff 

[' (kpc)] 


NGC 147 
NGC 185 
NGC 205 


00:33:12 
00:38:58 
00:40:22 


+48:30:31 
+48:20:14 
+41:41:07 


dE5 
dE3 
dE5 


675 ± 27 
616 ± 26 
824 ± 27 


9.52 ± 0.07 
9.18 + 0.05 
8.07 + 0.07 


0.57 
0.61 
0.21 


-15.5 
-15.7 
-16.7 


0.44 
0.23 
0.50 


2.03 (0.40) 
1.50 (0.27) 
2.17 (0.52) 



Note. — The right ascension, declination, and morphol ogical type are take from N ASA/IPAC Extragalactic 
Database (NED). We adopt the distances determined by McConnachic ct al. (2005). The total apparent V- 
band magnitude for each dE is taken from the RC3 catalog. The absolute magnitude is calculated assuming 
an extinctio n values from ISchlegel et alJ l|1998n. Th e ellipticity is the ave rage value between l' — 5' from 
IKentl 11987ft . The effective, half-light radii are from [de~Rijcke et all 1 )20061) , based on near-infrared 2MASS 
photometry. 



TABLE 2 

Keck/DEIMOS Multi-Slitmask Observing Parameters 



Mask 
Name 


a (J2000) 
(h : m : s) 


<5 (J2000) 
(•:':") 


PA 
(deg) 


^cxp 

(sec) 


# of slits 


% useful 
spectra 


N147.1 


0:32:54.2 


48 


22:48.0 


34 


3300 


190 


87% 


N147.2 


0:32:56.6 


48 


22:48.0 


31 


3600 


184 


79% 


N147.3 


0:33:35.0 


48 


37:09.5 


34 


3600 


195 


72% 


N147.4 


0:33:36.0 


48 


37:09.5 


214 


3600 


183 


82% 


N185.1 


0:39:36.0 


48 


26:06.0 


41 


3240 


155 


92% 


N185.2 


0:39:14.4 


48 


26:24.0 


221 


2040 


143 


78% 


N185.3 


0:38:39.5 


48 


14:58.2 


41 


3120 


135 


88% 


N185.4 


0:38:25.4 


48 


15:42.1 


221 


3060 


139 


42% 


N185.5 


0:39:03.6 


48 


26:24.0 


41 


3060 


143 


92% 


Note. - 


— Right ascension, declination, position angle and total 


exposure 



time for each Keck/DEIMOS slitmask. The final two columns refer to the 
total number of slitlets on each mask and the percentage of those slitlcts for 
which a rcdshift was measured. Note that the percentage of good slits for 
N185_4 is low due to a data reduction failure for 4 out of 8 CCDs. 
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TABLE 3 
Observed and Model Quantities 



Row 


Quantity 


Units 


NGC 147 


NGC 185 


Photometric Properties 






(1) 




mag 


s -I- n 
o.yo zll u.uo 


s ^7 -t- n 


(2) 


(B - V) 


mag 


0.76 ± 0.07 


0.72 ± 0.07 


(3) 


PA 


degrees 


208.4 


222.9 


(1) 


qa 




0.56 


0.77 


(->) 


i 


degrees 


77.5 


63.7 


«>) 






0.53 


0.70 


(7) 


n 




1.04 


1.76 


(«) 


jo(it/q a ) 


L Q pc- 3 


1.074 x 10~ 2 


3.008 x 10- 2 


(9) 


b 


arcmin 


10.72 


5.48 


(10) 


a 




-0.940 


-1.439 


(11) 


5 




-3.982 


-2.038 


(12) 


Ax 


arcmin 


1.87 


2.00 


(13) 


Ay 


arcmin 


5.00 


8.00 


Dynamical Properties 






(14) 


Vsys 


km s 


-193.1 ± 0.8 


-203.8 ± 1.1 


(15) 


^max,obs 


km s — 1 


17 ± 2 


15 ± 5 


(16) 




km s — 1 


16 ± 1 


24 ± 1 


(17) 


(Vrms) 


km s — 1 


19.1 ± 1.1 


24.6 ± 1.2 


(18) 


k 




0.74 ±0.06 


0.47 ± 0.08 


(19) 


M/L v 




4.24 ±0.56 


4.57 ± 0.56 


(20) 


M/L B 




4.71 ±0.69 


4.90 ± 0.68 


(21) 


Mass 


M 


5.6 x 10 8 


7.2 x 10 8 



(22) [Fc/H] dex -1.1 ±0.1 -1.3 ±0.1 

(23) cr [Pc/H] dex 0.4 0.5 

NOTE. — Quantities for the two dE galaxies, with uncertainties included where 
relevant. Flow (1) lists the observed total my Q magnitude from do Vaucoulours 
ct al. (1991; RC3), corrected for the foreground extinction listed in Table 1. 
Row (2) lists the extinction corrected color (B - V)q from RC3. Row (3) gives 
the adopted major-axis position angle, based on the outermost radius in Kent's 
(1987) photometry. Row (4) lists the apparent axial ratio of the models, which 
is the average axial ratio in the radial range 1 - 5' in the data of Kent (1987). 
Row (5) lists the statistically most likely inclination of each galaxy, based on 
the observed axial ratio q a and the distribution of intrinsic axial ratios from 
Tremblay Sz Merritt (1995). Row (6) lists the true axial ratio of each galaxy, 
based on the observed axial ratio q a and the listed inclination. Row (7) lists the 
value n of the best-fitting Sersic profile, which is shown in Figure^ Rows (8)- 
(11) give the parameters of the best-fitting luminosity density parameterization 
j(R, z) in equation Q. Row (12) lists the width Aa; of the bins used in binning 
the stellar velocities along the major-axis. Row (13) lists the approximate extent 
in the minor-axis direction of the region around the major-axis over which data 
was obtained. Row (14) lists the heliocentric systemic velocity V sys . Row (15) 
lists the weighted average value of V rms inside 4 7 . Row (16) lists the maximum 
observed rotational velocity, and row (17) lists the averaged observed velocity 
dispersion. Rows (18)-(20) arc the mass-to-light ratio M/L, and the rotation 
parameter k (see eq. \2^) inferred by fitting the models to our new kinematical 
data. The listed uncertainties in M j L include the propagated uncertainties in 
the galaxy distances D and the photometric uncertainties in the surface bright- 
ness profiles and color transformations. The S-band mass- to- light ratios were 
obtained from the modeled V-band values using the galaxy color (B - V)q and 
the solar color (B - V)q = 0.65 (Binncy &z Mcrrifield 1998). The characteristic 
mass M = (M/L) X L of each gala xy is listed in row (21). Row (22) lists the 
average metallicity as described in $ f3~2l row (23) lists the internal mctallicity 
dispersion in each dE. 
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TABLE 4 
Major-Axis Profile of NGC 147 



X 


a (J2000) 


S (J2000) 
(":':") 


V 


Verr 






[Fc/H] 


[Fc/H] orr 


(arcmin) 


(h : m : s) 


( km s _1 ) 


( km s — 1 ) 


( km s —1 ) 


(km s" 1 ) 


(dex) 


(dex) 


-14.7 


00:32:34.5 


+48:19:25.3 


-15.4 


3.3 


14.2 


2.4 


-0.7 


0.1 


-13.2 


00:32:38.6 


+48:20:35.5 


-18.4 


3.2 


16.4 


2.3 


-0.9 


0.1 


-11.1 


00:32:44.0 


+48:22:08.8 


-16.4 


3.2 


16.7 


2.3 


-0.7 


0.1 


-9.4 


00:32:48.6 


+48:23:27.6 


-11.0 


3.5 


21.5 


2.5 


-0.8 


0.1 


-7.5 


00:32:53.4 


+48:24:50.4 


-11.5 


3.4 


18.0 


2.5 


-0.9 


0.1 


—5.6 




4-4K-9R-1 7 9 


— 10.0 


2.7 


15.5 


2 


—0.7 


1 


-3.7 


00:33:03.5 


+48:27:44.3 


-6.8 


3.8 


21.3 


2.8 


-1.0 


0.1 


-1.9 


00:33:08.3 


+48:29:07.1 


-2.6 


3.6 


21.1 


2.7 


-1.1 


0.2 


-0.1 


00:33:12.8 


+48:30:25.9 


-1.8 


2.8 


17.2 


2.1 


-0.8 


0.2 


1.9 


00:33:18.1 


+48:31:57.0 


1.9 


3.7 


18.3 


2.9 


-1.0 


0.2 


3.7 


00:33:22.9 


+48:33:19.8 


8.7 


2.9 


16.2 


2.2 


-0.6 


0.1 


5.6 


00:33:27.8 


+48:34:44.8 


12.9 


3.7 


20.5 


2.7 


-1.1 


0.1 


7.4 


00:33:32.7 


+48:36:09.4 


11.8 


2.8 


15.6 


2.1 


-0.9 


0.1 


9.2 


00:33:37.4 


+48:37:30.4 


14.8 


2.8 


16.4 


2.1 


-0.9 


0.1 


11.3 


00:33:42.7 


+48:39:01.4 


13.8 


2.5 


10.1 


1.9 


-1.2 


0.2 


13.1 


00:33:47.5 


+48:40:24.2 


16.3 


2.8 


9.2 


2.2 


-1.3 


0.2 


14.7 


00:33:51.8 


+48:41:38.8 


17.2 


3.5 


7.7 


3.2 


-0.8 


0.1 


Note. — 


The velocity (V), velocity dispersion (a) 


and mctallicity ([Fc/H]) as 


a function of the binned 


major-axis 



distance x in arcminutcs away from the galaxy center. The right ascension and declination are given for each radial 
bin: positive radial bins on the North-East side of the galaxy, negative radial bins are on the South- West. 



TABLE 5 
Major-Axis Profile of NGC 185 



X 

(arcmin) 


a (J2000) 
(h : m : s) 


<5 (J2000) 
(":':«) 


V 

(km s" 1 ) 


Verr 

(km s" 1 ) 


(7 

(km s" 1 ) 


Cerr 

(km s" 1 ) 


[Fe/H] 
(dex) 


[Fc/H] crr 
(dex) 


-12.3 


00:38:22.2 


+48:11:53.5 


-12.6 


9.1 


28.8 


6.5 


-1.4 


0.1 


-9.9 


00:38:29.0 


+48:13:30.0 


-9.3 


5.1 


24.0 


3.9 


-1.4 


0.1 


-8.1 


00:38:34.5 


+48:14:46.0 


-7.0 


3.6 


19.6 


2.6 


-0.8 


0.1 


-5.9 


00:38:40.6 


+48:16:10.9 


-3.6 


6.7 


27.6 


1.8 


-1.2 


0.2 


-4.0 


00:38:46.2 


+48:17:29.4 


-7.9 


6.8 


31.8 


1.9 


-0.9 


0.3 


-2.0 


00:38:52.1 


+48:18:52.2 


-6.1 


4.5 


26.5 


3.2 


-0.7 


0.3 


-0.1 


00:38:57.8 


+48:20:12.8 


0.5 


3.8 


23.9 


2.7 


-0.8 


0.2 


1.9 


00:39:03.6 


+48:21:33.5 


-0.3 


3.1 


23.0 


2.2 


-0.9 


0.1 


1.0 


00:39:09.7 


+48:22:58.4 


3.3 


3.3 


23.6 


2.4 


-0.9 


0.1 


6.0 


00:39:15.6 


+48:24:21.2 


7.1 


3.2 


20.5 


2.3 


-1.4 


0.1 


8.0 


00:39:21.4 


+48:25:41.9 


5.8 


3.3 


20.9 


2.4 


-1.3 


0.1 


10.0 


00:39:27.1 


+48:27:02.2 


11.7 


3.9 


20.0 


2.8 


-1.4 


0.1 


11.9 


00:39:32.9 


+48:28:22.8 


11.6 


4.7 


18.8 


3.5 


-1.5 


0.1 


14.4 


00:39:40.2 


+48:30:05.8 


10.5 


6.1 


22.0 


4.1 


-1.5 


0.1 



Note. — See comments for Table 4. 



TABLE 6 

Keck/DEIMOS Velocity Measurements for NGC 147 



i Name a (J2000) <5 (J2000) /-mag (R - I) v v cr r 

(li m s) (° ' ") ( ma g) ( ma g) (km s _1 ) (kms -1 ) 



1 


11231 


00:33:02.4 


+48:23:16.7 


21.3 


0.97 


-235.5 


3.3 


2 


4456 


00:32:23.9 


+48:19:53.8 


20.5 


1.09 


-235.5 


3.3 


3 


4608 


00:32:25.8 


+48:18:57.2 


20.7 


1.14 


-206.2 


2.8 


4 


4778 


00:32:27.6 


+48:18:13.8 


20.9 


0.76 


-203.4 


1.9 


5 


4894 


00:32:29.7 


+48:19:03.8 


21.1 


0.94 


-209.9 


3.2 



Note. — Velocity measurements for member stars of NGC 147. Position, apparent 
/—band magnitude, {R — I) color, heliocentric radial velocity (v), and velocity error (v eTT ) 
for each star as determined in § |2.3I This table will be published in its entirety in the 
electronic edition of the Astrophysical Journal. 
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TABLE 7 

Keck/DEIMOS Velocity Measurements for NGC 185 



i Name a (J2000) <5 (J2000) /-mag (R - 1) v v CII 

(h m s) (° ' ") (mag) (mag) (kms" 1 ) (kms" 1 ) 



1 


11396 


00:38:58.6 


+48:21:52.3 


20.5 


0.66 


-203.3 


3.8 


2 


11487 


00:38:58.9 


+48:23:33.3 


20.9 


1.18 


-202.3 


4.4 


3 


11515 


00:38:59.0 


+48:22:03.4 


20.8 


1.01 


-190.0 


5.9 


1 


11522 


00:38:59.0 


+48:22:34.5 


21.0 


0.69 


-215.6 


4.1 


5 


12017 


00:39:00.2 


+48:22:09.0 


20.1 


0.96 


-175.7 


2.8 



Note. — 



Same as Table 6, but for NGC 185. 



